Non-Gaussian Radio- Wave Scattering in the Interstellar Medium 
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ABSTRACT 

It was recently suggested by Boldyrev & Gwinn that the characteristics of radio scintillations 
from distant pulsars are best understood if the interstellar electron-density fluctuations that cause 
the time broadening of the radio pulses obey non-Gaussian statistics. In this picture the density 
fluctuations are inferred to be strong on very small scales (~ 10 8 — 10 10 cm). We argue that such 
density structures could correspond to the ionized boundaries of molecular regions (clouds) and 
demonstrate that the power-law distribution of scattering angles that is required to match the 
observations arises naturally from the expected intersections of our line of sight with randomly 
distributed, thin, approximately spherical ionized shells of this type. We show that the observed 
change in the time-broadening behavior for pulsar dispersion measures < 30pccm -3 is consistent 
with the expected effect of the general ISM turbulence, which should dominate the scattering 
for nearby pulsars. We also point out that if the clouds are ionized by nearby stars, then their 
boundaries may become turbulent on account of an ionization front instability. This turbulence 
could be an alternative cause of the inferred density structures. An additional effect that might 
contribute to the strength of the small-scale fluctuations in this case is the expected flattening of 
the turbulent density spectrum when the eddy sizes approach the proton gyroscale. 

Subject headings: ISM:general-pulsars:general-scattering-MHD-turbulence 



1. Introduction 

Radio signals received from distant pulsars fluc- 
tuate in time and space due to inhomogeneities in 
the ionized component of the interstellar medium. 
Observations of pulsar scintillations have thus 
served as a valuable tool for reconstructing the 
statistics of interstellar electron-density fluctua- 
tions (e.g., Rickett 1990; Armstrong et al. 1995; 
Scalo & Elmegreen 2004) . Among a variety of ob- 
servational quantities, special attention has been 
given to the time shape of the pulsar signal I(t) 
(the signal intensity as a function of time). This 
function arises from the interference of waves prop- 
agating along different paths (similar to the Feyn- 
man interpretation of the propagation of a quan- 
tum particle). The deviation of the wave paths 
from a straight line causes a time-broadening of 
the arriving signal. The intensity function I(r) 
can be interpreted as the probability density func- 



tion (PDF) of the time delays caused by different 
paths of propagation. Since the statistics of the 
path deviations are directly related to the statis- 
tics of electron density fluctuations, the function 
J(r) provides information on the latter. 

Boldyrev & Gwinn (2003a,b, 2005) recently 
used the observed time shapes of the measured 
pulses to probe the shape of the distribution func- 
tion of the underlying density fluctuations. Based 
on a comparison of analytical results and observa- 
tional data, they proposed that the density fluc- 
tuations responsible for the pulse broadening are 
non-Gaussian, with the distribution function of 
the resulting angular path deviations having a 
slowly declining, power-law asymptotic. They fur- 
ther deduced that the density fluctuations should 
be strong on very small scales (estimated in § 3 to 
be ~ 10 s - 10 10 cm). 

After being suitably rescaled by the pulsar dis- 
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tance and the observation wavelength, the ob- 
served pulse time shapes exhibit a somewhat "uni- 
versal" behavior (e.g., Bhat et al. 2004; Lohmer 
et al. 2004), which indicates that similar pro- 
cesses govern the density fluctuations along differ- 
ent lines of sight. Two intriguing questions then 
arise regarding the physics of interstellar plasma 
fluctuations. First, how are strong small-scale 
density fluctuations generated? And second, how 
can density fluctuations produce a non-Gaussian, 
power-law distribution of scattering angles? 

In addressing these questions one can consider 
two distinct (although not mutually exclusive) 
possibilities. In one interpretation the observed 
universality is attributed to a general statisti- 
cal property of density fluctuations. For exam- 
ple, in Kolmogorov-type turbulence, small-scale 
fluctuations arise as a result of an energy cas- 
cade from large to small scales. The proper- 
ties of the observed signals might then be linked 
to the universal properties of such a cascade. 
The second possibility is that density fluctuations 
are strongly nonuniform and spatially intermittent 
(e.g., clumps, filaments, and shocks). In this case, 
the universality could reflect a certain inherent 
property of the density structures. This property 
should be fairly robust in view of the fact that dif- 
ferent mechanisms (not necessarily related to tur- 
bulence) could in principle give rise to a particular 
structure. For example, scattering by shocks may 
be universal if the only property that matters is a 
sharp density discontinuity (see, e.g., Boldyrev & 
Gwinn 2005) . 

Let us briefly discuss the first possibility. We 
assume that interstellar scattering occurs in tur- 
bulent ionized regions (HII regions) of mean den- 
sity n,2 ~ 10 2 cm~ 3 . Such regions appear when 
bright stars turn on within or in the vicinity of cold 
molecular clouds (e.g., Dyson & Williams 1997). 
Large-scale density fluctuations in a cold molec- 
ular cloud of mean density no ~ 10 2 cm -3 are 
of the order Ano ~ 10 2 cm~ 3 , since such clouds 
are usually turbulent with Mach numbers greater 
than 1 (Larson 1981). When cold neutral gas in 
a molecular cloud is ionized by a newborn mas- 
sive star emitting ~ 10 49 ionizing photons s _1 , 
the radius of the ionized sphere (the Stromgren 
radius) is R s ~ 10 19 cm. This ionization hap- 
pens very fast, during the first 10 4 years after the 
radiation turns on. This time interval is much 



shorter than the sound crossing time in the ion- 
ized gas within the Stromgren sphere, where the 
sound speed is Ci ~ 10 6 cm s" 1 . Subsequently, 
the initial density inhomogencitics relax due to 
pressure gradients, which drive turbulent internal 
motions in the HII region. The longest relaxation 
time corresponds to the largest spatial scale, and 
therefore the turbulence that survives over the life- 
time (~ 10 6 yr) of the star is effectively stirred on 
the scale Z ut ~ Rs, where the velocity and den- 
sity fluctuations are of the order of Ci and An , 
respectively. 

The interstellar medium is magnetized and the 
magnetic field is likely to play a role in any tur- 
bulent energy cascade. As discussed by several 
authors (e.g., Higdon 1984; Goldreich & Sridhar 
1997; Lithwick & Goldreich 2001), density per- 
turbations in an MHD turbulence are associated 
with compressible and entropy modes and are pas- 
sively advected by the Alfvenic cascade toward the 
smallest (dissipative) spatial scales. The cascade 
proceeds predominantly in the direction perpen- 
dicular to the magnetic field and continues until 
the fluctuations reach the scale of the ion gyrora- 
dius pi. A typical strength of the magnetic field in 
the Galaxy is a few /jG (Zweibel & Heiles 1997). 
For an ion temperature T, ~ 8000 K and mag- 
netic field amplitude B ~ 3 pG, the gyroradius is 
Pi ~ 3 x 10 7 cm, so the density fluctuations can 
reach very small scales. 

However, a problem arises when we address the 
statistics of density fluctuations that are advected 
as a passive scalar in MHD turbulence. Numerical 
results show that passive scalar fluctuations in an 
incompressible turbulence have an exponential (or 
stretched exponential) distribution (e.g., Warhaft 
2000), not the power-law distribution that we are 
trying to explain. Another problem arises when 
we recall that the fluctuations were inferred to be 
strong on scales as small as ~ 10 8 — 10 10 cm. To 
reach these scales, the envisioned cascade would 
have to pass through a scale where the eddy 
turnover time is comparable to the radiative cool- 
ing time of the plasma, at which point the den- 
sity fluctuation may be strongly damped (Lith- 
wick & Goldreich 2001). The radiative cooling 
time for ionized gas of temperature ~ 8000 K is 
given by t coo i ~ 6 x lO 11 ?^ 1 s (Spitzer 1978), so for 
77,2 ~ 10 2 cm -3 and a Kolmogorov-likc turbulence, 
the cooling scale is i coo i ~ 10 14 cm, significantly 
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larger than the scales we are trying to account 
for. 

In view of the foregoing arguments it appears 
difficult to attribute the inferred non-Gaussian 
scattering by small-scale density fluctuations to a 
nearly incompressible MHD turbulence excited on 
scales Z ou t ~ 10 pc. In § 5 we discuss a possible 
means of exciting the turbulence on scales smaller 
than Z coo i (to account for the inferred strength of 
the small-scale fluctuations) and appeal to kinetic 
effects (which could make the turbulence com- 
pressible) for further enhancement of the fluctu- 
ations and the possible establishment of a power- 
law spectrum on scales close to pj. However, our 
main focus in this paper is on exploring the sec- 
ond of the above two possibilities for the origin 
of the inferred density fluctuations, namely, that 
the spatial distribution of the scattering electrons 
is highly intermittent. Specifically, we propose 
that these electrons are contained in thin, shell- 
like structures. Although we show that the sta- 
tistical properties of the scattered radio pulses are 
independent of the physical origin of the shells, 
we argue that the latter can be most plausibly 
identified with the thin surface layers of molecular 
regions (clouds) that become ionized by the inter- 
stellar radiation field or by ionizing photons from 
nearby stars. 

The ionized surfaces of dense molecular clouds 
were previously invoked by Yusef-Zadeh et al. 
(1994) to explain (on the assumption that these 
layers are magnetized and turbulent) the observed 
anisotropy in the angular broadening of the radio 
source Sgr A* in the Galactic center. Furthermore, 
localized scattering structures (such as shocks or 
filaments) have been invoked by Lambert & Rick- 
ctt (2000) and by Cordes & Lazio (2001) to explain 
the observed scaling of pulsar signals with fre- 
quency and the results of certain diffractive mea- 
surements. Small-scale localized electron density 
structures have also been inferred from observa- 
tions of extreme scattering events (e.g., Lazio, Fey 
& Gaume 2001), pulsar arclet structures (e.g., Hill 
et al 2005), and pulsar scintillations (e.g., Cordes, 
Weisberg, & Boriakoff 1985). Although these in- 
vestigations were carried out in the context of the 
Gaussian model and thus do not bear directly on 
the results presented here (which aim to justify 
a non-Gaussian interpretation), their appeal to 
sharply delineated scattering regions has inspired 



our work. 

In this paper we demonstrate that the non- 
Gaussian, power-law statistics inferred for the an- 
gular scattering of pulsar radio pulses need not be 
related to the character of density fluctuations in- 
side turbulent scattering regions but can be fully 
explained in terms of the spatially intermittent 
nature of thin ionized shells in the interstellar 
medium (which in principle need not even be tur- 
bulent). Using purely geometric-optics considera- 
tions, we show that scattering by thin spherical 
shells produces a distribution of angular devia- 
tions A9 with an asymptotic power-law probabil- 
ity density P(A9) cx A0~ 5 / 3 . Interestingly, this is 
precisely the form deduced in Boldyrev & Gwinn 
(2003a) by applying the non-Gaussian (Levy) scat- 
tering theory to the pulsar scintillation data. This 
correspondence strongly suggests that the scatter- 
ing regions that produce refractive radio scintilla- 
tions have a morphology of curved shells. 

After briefly outlining the theory of pulsar scin- 
tillations (§ 2), we discuss the key differences be- 
tween Gaussian and non-Gaussian models of den- 
sity fluctuations (§3). In § 4 we flesh out our pro- 
posal that scattering by thin ionized boundaries 
of molecular clouds can naturally explain the in- 
ferred non-Gaussian statistics of the scintillations. 
In § 5 we consider the case in which the ionizing 
radiation is strong enough to trigger an ionization 
front instability and possibly the development of 
turbulence within the ionized cloud boundaries. 
Our results are summarized in § 6. 

2. Pulsar Scintillations and Interstellar 
Density Fluctuations 

In this section we explain how the statistics of 
electron-density fluctuations can be diagnosed by 
observations of the time shapes of pulsar signals. 
To this end, we assume that the scattering (refrac- 
tion) of a radio wave coming from a distant pul- 
sar occurs in separate planar phase screens, uni- 
formly placed along the line of sight to the pul- 
sar. We denote the line-of-sight axis by z and the 
coordinates in the perpendicular plane by a two- 
dimensional vector y. We also assume that the 
scattered wave is planar and propagates close to 
the line of sight. To analyze this model one can 
make use of the so-called parabolic approximation, 
which assumes that the wave amplitude changes 



3 



slowly compared to its phase. The details of this 
analysis can be found in a number of references 
(e.g., Tatarskii 1961; Uscinski 1974; Williamson 
1972; Lee & Jokipii 1975; Rickett 1990; Boldyrev 
& Gwinn 2005); here we only give a qualitative 
discussion of the results. 

Each scattering screen generates phase fluctua- 
tions that cause the transmitted wave to contain 
many different angular components in its spec- 
trum rather than follow a single direction. The 
spatial scale of the fluctuations that contribute to 
an angular deviation AO (expressed in radians) is 
y <~ X/A8, where A is the wavelength. The dif- 
ference of phases accumulated in a single screen 
along two lines of sight separated by the distance 
y is given by A$ = Ar AA/"(y). It is proportional 
to the line-of-sight integrated density difference 



AA%) 



Kyi, 2) - n{y 2 ,z))dz , 



(1) 



where y = yi — y 2 and r = e 2 /(m e c 2 ) — 2.82 X 
10~ 13 cm is the classical electron radius. The inte- 
gration length l should exceed the characteristic 
size of the scattering region. To produce an angu- 
lar deviation A9, the phase difference should be 
A$ ~ 2ir on the scale y <~ X/A8. (where here, as 
in the rest of this paper, angles are measured in 
radians) . 

The signal emitted by a pulsar has a rather 
short duration (several percent of the pulsar pe- 
riod). The received signal is broader, with its in- 
tensity I(t) exhibiting a sharp rise and a slow de- 
cline. This happens because waves propagating 
along different paths from the pulsar to Earth ex- 
perience different angular deviations by the elec- 
tron density prisms that intercept them. This, 
in turn, leads to different time delays relative to 
straight (undeflected) propagation. The signal in- 
tensity corresponding to a time delay r is propor- 
tional to the probability for a path to be delayed 
by t (e.g., Boldyrev & Gwinn 2005). 

Assume that the number of scattering screens 
is N and the distance between two neighboring 
screens along the line of sight is z , so that d = 
Nzo is the distance to the pulsar. The angle of the 
path in the segment between screens m — 1 and m 
is a sum of the angular increments accumulated 
at the preceding segments, m = J2T=i ^ n 
general, A0 is a two-dimensional vector, pointing 



in the direction of the path deviation. The time 
delay introduced by the mth scattering segment is 
therefore Ar m — (z /c)(l — cos6* m ) ~ (z o /2c)0. 
and the total time delay is the sum of the individ- 
ual delays, 
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N 
m—l 



(2) 



In this equation the individual angular increments 
A0 are independently and identically distributed 
two-dimensional vectors. The intensity I(t) of the 
received signal, which is effectively the PDF of r, 
can be found numerically using equation (2) once 
the distribution function for the vectors A0 S is 
known. 

In the classical theory of scintillations the vec- 
tors A0 S are assumed to have an independent 
and identical Gaussian distribution, in which case 
an exact analytic solution for the r-distribution 
function I(t) can be derived (e.g., Uscinski 1974; 
Williamson 1972). However, the scaling and 
shapes of pulse profiles predicted by this theory 
seem to disagree with the observational results for 
distant pulsars, as first noted by Sutton (1971) 
and Williamson (1974). In a recently proposed 
alternative model, Boldyrev & Gwinn (2003a) 
demonstrated that the observed profiles can be 
matched if the individual fluctuations A6 S have a 
non-Gaussian, power-law-like declining distribu- 
tion. In this case the wave path angle 9 m exhibits 
a so-called Levy flight rather than the standard 
random walk. 

The Levy distribution function is a general dis- 
tribution function whose convolution with itself 
produces the same function again (appropriately 
rescaled). The Gaussian distribution is a par- 
ticular case of a Levy distribution. A detailed 
discussion of Levy distributions can be found in 
Klafter et al. (1995); the application to wave 
propagation in the interstellar medium is consid- 
ered in Boldyrev & Gwinn (2003a,b, 2005). For 
our present purposes, the main relevant property 
of the Levy distribution is its slow (power-law) 
asymptotic decline. Whereas 9 m oc m for a Gaus- 
sian random walk, the scaling produced by a Levy 
flight is 9 2 n cx m 2 /' 3 , where (3 (which satisfies 
< (3 < 2) is the parameter of the Levy distribu- 
tion. For (3 < 2, the Levy PDF exhibits a power- 
law decline, Pp(A0) cx lA^- 1 -^ for |A0| > A0 O , 
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where A6> is a typical angular fluctuation value. 

Pulsar signals are observed to satisfy the scal- 
ing r oc TV 4 (more precisely, the observationally 
inferred scaling is r oc DM 4 , where the pul- 
sar dispersion measure, DM, is proportional to 
the distance to the pulsar and, therefore, to the 
number of scattering events, N), which motivated 
Boldyrev & Gwinn (2003a, b, 2005) to propose that 
radio scintillations are described by a Levy model 
with (5 ~ 2/3. A Gaussian distribution corre- 
sponds to (3 = 2 and implies a linear scaling of 
Q 2 ^ with m, which is not consistent with the ob- 
servational data. In § 3 we demonstrate that the 
difference in the implied scaling between these two 
models corresponds to a profound difference in the 
physics of the respective scattering mechanisms. 

3. Levy Statistics of Scintillations 

We assume that the PDF of AO has a charac- 
teristic width A6* (corresponding, say, to the 1/e 
amplitude level), and we denote the width of the 
resulting PDF of t by r . Equation (2) and the 
scaling Q 2 m oc m imply that r scales as (A6 ) 2 N 2 
in a Gaussian model. On the other hand, if we 
consider To to arise from the action of individual 
fluctuations, each of magnitude AO, we find (us- 
ing eq. [2] again) that To oc (A6*) 2 iV 3 . By compar- 
ing these two expressions we obtain, as expected, 
A9 ~ AOq/N 1 / 2 , which is consistent with the no- 
tion that the typical time delay may be produced 
by many small deflections. 

The situation is quite different in the Levy 
case with (3 — 2/3 < 1. Equation (2) implies 
t oc (A8 ) 2 N 4 , and if we compare this expression 
with our estimate in terms of individual fluctua- 
tions, To oc (A8) 2 N 3 , we infer that such fluctu- 
ations must be large, A9 ~ N 1 / 2 A9o- But such 
large individual fluctuations have a small proba- 
bility, P oc l/{AB) 1+fi oc 1/N 5 / 6 , and it is there- 
fore highly unlikely that all N individual fluctua- 
tions are of that order. Instead, the most proba- 
ble situation is that one of the fluctuations is ex- 
tremely large, AO = A6 max ^ N 3 / 2 A9 , whereas 
all the others are small (< A9o). Thus, in each 
particular realization, the total scattering is most 
probably produced by only one, randomly chosen, 
screen. This behavior is a general property of sums 
of Levy-distributed variables with (3 < 1 (Feller 
1971). When the pulse shape is averaged over a 



certain time interval (or over a statistical ensem- 
ble), these screens can be different in different re- 
alizations, since the interstellar medium fluctuates 
and both the pulsar and the Earth move through 
it. 

We can now estimate the width of the pulse 
as t ~ z iV(A6» max ) 2 /c ~ d(A8 max ) 2 /c, where 
d = Nz is again the pulsar distance. This ex- 
pression has a geometric meaning, namely, that 
only one refractive event dominates the path devi- 
ation from a straight line. For the typical param- 
eters of distant pulsars (e.g., Bhat et al. 2004), 
d ~ 3 x 10 22 cm, t ~ 0.01s, and A ~ 30cm 
(corresponding to a frequency of 1GHz), we ob- 
tain A9 

max r " J 

10 7 rad and a density fluctuation 
scale y ~ X/ A9 max ~3x 10 8 cm. In view of the 
scaling r a d 4 , closer pulsars (characterized by 
d <~ 3 x 10 21 cm and t ~ 10~ 6 s) correspond to 
A8 max - 3 x 10" 9 rad and y ~ 10 10 cm. 

What regions in the interstellar medium might 
be responsible for the indicated strong scattering 
events? To address this question, we first consider 
whether the maximal angular deviation estimated 
above could be produced by regions of uniform and 
homogeneous turbulence. In the following discus- 
sion we therefore neglect the contribution of the 
boundaries of the turbulent regions or of any other 
structures. 

The amplitude of the density fluctuations that 
give rise to the large inferred angular deviations 
can be estimated from the results presented after 
equation (1), which imply 

AAf(y)/y ~ A6 max /{\ 2 r ) ~ 10 3 cm" 3 . (3) 

This defines the value that the density fluctuations 
must be able to reach on the scale y. Suppose that 
the density fluctuations are associated with a tur- 
bulent cascade. If the density is passively advected 
by the fluid turbulence (e.g., Lithwick & Goldre- 
ich 2001) then the density scaling follows that 
of the velocity field, and for a spatially homoge- 
neous turbulence with an outer scale / ut 

we have 

Sn(l)/An ~ Sv{l)/Av ~ (Z/Ut) Q/2 . Here Av 
and Ano are, respectively, the rms values of the 
velocity and density fluctuations on the scale l out 
and a is the turbulence scaling exponent, which 
ranges between 2/3 (the Kolmogorov case) and 2. 
This estimate neglects damping effects, which can 
be different for 5n and 8v; we discuss damping 
further on in this section. 
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As we pointed out in § 1, the distribution 
function of density fluctuation amplitudes in Kol- 
mogorov turbulence declines faster than a power 
law. The central limit theorem then implies that 
the distribution of line-of-sight integrated density 
fluctuations is Gaussian. Under these conditions, 
the y-dependence of the integrated density differ- 
ence is 

AA%) ~ zi /2 ^ t 2 An (y/Ut) 5 , (4) 

where 6 = (l + a)/2 for a < 1 and 6 = 1 for a > 1. 
Equation (4) follows directly from the expression 
for the correlator ([AW(y)] 2 ) in a homogeneous 
and isotropic turbulence. Note that this correlator 
is proportional to the integration distance lo, as 
expected for Gaussian fluctuations. 

In the standard picture of scintillations (the 
first scenario discussed in § 1), turbulent HII re- 
gions around bright stars are considered to be the 
likely sites of radio- wave scattering. In this case 
the outer scale of the turbulence, l out ~ 10 19 cm, is 
of the order of the size lo of the Stromgrcn sphere. 
Setting Iq ~ l out in equation (4) and comparing 
with equation (3), we obtain the following condi- 
tion on the large-scale density fluctuations: 

An ~10 3 (y/Ut) (1 - a)/2 cm- 3 . (5) 

Taking the velocity difference at Z out to be of the 
order of the sound speed in the ionized region 
(C2 ~ 10 6 cm/s) and assuming (conservatively) 
that the turbulence in the H II region indeed 
obeys a Kolmogorov scaling (a = 2/3) and that 
y ~ 10 s cm, we find that equation (5) is satis- 
fied on this scale for electron density fluctuations 
An > 30 cm- 3 . 

Although this value is consistent with the typ- 
ical densities of molecular clouds, in reality the 
required value of Ano should be much larger (and 
therefore implausible) in view of the fact that the 
compressible and entropy MHD modes associated 
with the density fluctuations arc damped when 
the wave cascade passes through the cooling scale 
(see Lithwick & Goldreich 2001), which is about 
10 14 cm in this case (see § 1). 

The cooling length constraint limits the abil- 
ity of the standard picture to explain the strong, 
small-scale density fluctuations that are required 
in the Levy-flight scenario. Furthermore, as dis- 
cussed in § 1, the distribution of passive scalar 



fluctuations in the incompressible MHD model is 
not power-law-like. These two difficulties suggest 
that non-Gaussian scattering cannot be explained 
within the standard framework. In § 4 we demon- 
strate that such scattering may be explained if it is 
produced in thin shells — the ionized boundaries 
of molecular regions (clouds). 

4. Non-Gaussian Scattering by Ionized 
Boundaries of Molecular Clouds 

A typical value for the photon flux of ionizing 
interstellar radiation (hv > 13.6 cV) through a 
unit surface can be estimated by converting the 
standard Habing flux (Habing 1968) of dissociat- 
ing photons (6 < hv < 13.6 eV) using a black- 
body spectrum of temperature T~3x 10 4 K (see 
Yusef-Zadeh et al. 1994). This gives J ~ 2x 
10 7 cm~ 2 s _1 . In the vicinity of a bright star the 
ionizing flux can be several orders of magnitude 
higher. For example, at a distance of 1 pc from a 
star emitting S = 10 49 ionizing photons this 
flux is J ~ 10 11 cm _2 s _1 . The mean free path 
of an ionizing photon that penetrates a neutral 
medium of density no is A p h ~ l/(aono), where 
ao = 6.8 x 10~ 18 cm 2 is the hydrogen ionization 
cross section at the threshold energy hv = 13.6 eV 
(e.g., Spitzer 1978). 

When a homogeneous molecular cloud is irra- 
diated by an ionizing photon flux J it develops an 
ionized "skin" whose width Ar^ can be found from 
the ionization balance condition J = fan^/S.ri, 
where fa is the recombination coefficient (fa ~ 
2 x 10~ 13 cm 3 s- 1 for T ~ 10 4 K; e.g., Dyson & 
Williams 1997). For n ~ 10 2 cm -3 and J = J 
one finds Ar^ ~ 10 16 cm, which can be taken as the 
thickness of the ionized boundary. The width of 
the transition layer between the ionized skin and 
the neutral interior of the cloud is of the order of 
Aph (~ 10 15 cm for the adopted value of no). 

We now demonstrate that thin scattering shells 
(representing the ionized boundaries of molecular 
clouds) can produce non-Gaussian scintillations 
even if they are not turbulent. In particular, we 
show that scattering by thin curved layers (which 
we represent for simplicity as thin spherical shells) 
produces a power-law distribution of the scatter- 
ing angles. Most remarkably, the derived power- 
law index coincides with the value originally in- 
ferred by Boldyrev & Gwinn (2003a) from an in- 
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Fig. 1. — Scattering by a thin spherical shell. The 
ray's angular deviation is very small, so its trajec- 
tory is well approximated by a straight line. 

terpretation of the scintillation data in terms of 
the Levy theory (see § 2 and 3). The agreement 
with the observations suggests that a shell-like 
structure is a general feature of the ionized scat- 
tering regions in the interstellar medium. 

To derive these results, we assume that each 
of the electron-scattering regions has the form 
of a thin spherical shell of radius R and thick- 
ness AR <C R (see Fig. 1). We denote by An the 
electron-density excess in the shell relative to its 
surroundings and employ geometric optics to cal- 
culate the angular deviation A9 experienced by a 
ray that intercepts the shell at a distance b from 
the center. For a typical density excess An ~ 
10 2 cm -3 this deviation is very small, A9 <C 1, so 
the trajectory of the ray through the scattering 
region can be approximated by a straight line. 

Scattering by a shell may be viewed as a super- 
position of scatterings by two spheres of radii R 
and R—AR, respectively. The total angular devia- 
tion produced by a sphere of radius r is readily ob- 
tained using Snell's law and is given as a function 
of the impact parameter b by 9 r (b) = A9 s (b/r)(l — 
b 2 /r 2 ) -1 / 2 . In this expression A9 S = A 2 r An/7r, 
where A and r are again the wavelength and the 
classical electron radius, respectively. 

The angular deviation produced by the shell 
is A6(b) = 6 R (b) - 6 R - AR (b). It is given, up 
to first order in the small parameter AR/R, by 




AG AG C AG 

Fig. 2. — A sketch of the angular-deviation PDF 
on a log-log scale. 

A0(b) = A9 s (bAR/R 2 )(l-b 2 /R 2 )- 3 / 2 . Large an- 
gular deviations, A9(b) ^> A9 S AR/ R, correspond 
to b « R, and in this limit we obtain 

A8(b) » A9 s (AR/R)(2Ab/R)-^ 2 , (6) 

where Ab = R — b. When b is not close to R, the 
angular deviations are small, A9 ~ A9 S (AR/ R). 

The probability density for the ray to pass at 
a distance b from the center is estimated using 
P(b)db = 2irbdb/(irR 2 ) - (2/R)db, where we 
again assume b w R. The corresponding PDF of 
the deviation angle A9 can be found using equa- 
tion (6), which gives 

P(A9) ~ (2/3)(A^Ai?/i?) 2 / 3 (A0)- 5 / 3 . (7) 

This expression describes the asymptotic behavior 
of the PDF for A9 > A9 S AR/R. The derived ex- 
ponent (—5/3) of the power-law decline of P{A9) 
coincides with the exponent predicted in the Levy 
model (see § 2). Overall, the PDF P(A9) has a 
typical width of order A9 ~ A9 S (AR/R) (corre- 
sponding to impact parameters b not close to R). 
At larger angular deviations the PDF develops the 
power-law shape obtained above, which extends 
up to the cutoff value A9 C ~ A9 S (R/AR) 1 / 2 (at- 
tained when Ab decreases below ~ A_R). For even 
larger angular deviations (i.e., for Ab <C AR) the 
scattering is produced only by the outer sphere 
in Figure 1 and the PDF exhibits a fast decline, 
P(A9) oc (A9)- 3 . This PDF is sketched in Fig- 
ure 2. 

For our fiducial parameter values A9 S ~ 
lCT 8 rad and A9 max ~ l(T 7 rad (see § 3). 
Based on equation (6), the inferred value of 
the maximal angular deviation is attained when 
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(AR/R)(2Ab/R)^ 2 ~ A9 max /A6 s ~ 10. Given 
that Ab cannot be smaller than the shell thick- 
ness (i.e., Ab > AR), the shell thickness must 
satisfy AR/R < 1CP 2 . If the shell's radius of cur- 
vature is R ~ 10 pc, we have AR < 3 x 10 16 cm, 
which is in excellent agreement with our esti- 
mate of Ari in typical molecular clouds. We 
note, however, that to calculate the total prob- 
ability distribution function of the angular devi- 
ation AO one needs to convolve the distribution 
P(A8) (which was obtained for a boundary of 
given curvature and density) with the curvature 
radii and density distributions of real clouds. This 
could in principle change the value of the coeffi- 
cient (A9 S AR/R) 2 / 3 in equation (7) as well as the 
range of angular deviations over which the derived 
power-law asymptotic applies. The curvature radii 
of a cloud's boundary are typically smaller than 
the cloud size, so we may use the cloud size as an 
upper bound on the curvature R. Now, A8 S AR/R 
scales approximately as l/(noR). To the extent 
that the column density of molecular clouds is ap- 
proximately a constant (e.g., Larson 1981; McKee 
1989), the coefficient in equation (7) would change 
little from cloud to cloud. In a similar vein, the 
cutoff angular deviation A8 C scales approximately 
as i? -3 / 2 . Therefore, our choice of fiducial cloud 
radius (i? = 10 pc), which lies near the upper end 
of the observed distribution (e.g., Spitzer 1978), 
is conservative in that it leads to a lower bound 
on A9 C . All in all, we expect that integrating 
the expression (7) over a cloud distribution would 
have little effect on the derived asymptotic form 
of P(A8) and, in particular, would not change the 
inferred value of its power-law exponent. 

It is interesting to consider some observational 
implications of the proposed model. As we estab- 
lished in this section, AR/R < (A6 S / A8 max ) 2 . 
On the other hand, from § 3 we know that 
when the scattering angles have a Levy dis- 
tribution with (3 = 2/3, the number of scat- 
tering structures along the line of sight scales 
with the maximal deflection angle as iV 3 / 2 ~ 
A8 max /A8 Q ~ (A6 max /A8 S ){R/AR). There- 
fore, our considerations suggest that the number 
of scattering structures out to the characteristic 
distance d ~ 3 x 10 22 cm can be estimated as 
N > (A8 max / A0 S ) 2 ~ 100, or about one scatter- 
ing structure per 100 pc. This estimate is consis- 
tent with the overlap radius of radiatively cooled 



supernova shells (within which dense clouds are 
envisioned to form) in the McKee & Ostrikcr 
(1977) model of the interstellar medium. The 
estimated number of scattering regions could in 
principle be smaller (resulting in a weaker con- 
straint on their separation) since several indepen- 
dent scattering structures (ionized boundaries) 
might in practice form inside a single scattering 
region. 

The physical basis for the statistical proper- 
tics of the thin, curved scattering structures in 
the proposed interpretation is their intermittent 
spatial distribution. Specifically, the free elec- 
trons responsible for the scattering are envisioned 
to occupy only a small fraction of the interstel- 
lar medium — namely, the ionized boundaries 
of molecular clouds. Furthermore, the molecu- 
lar clouds themselves are sparsely distributed in 
space, so the effects of two different clouds on 
wave propagation are uncorrelated. Thus, a radio 
signal from a distant pulsar does not experience 
significant scattering over most of its propagation 
distance to Earth. However, when its trajectory 
intersects a curved, ionized cloud boundary, the 
scattering may be strong. As discussed in § 3, 
this behavior is characteristic of Levy-flight wave 
paths. The spatial intermittency in this scenario is 
distinct from the intermittent nature of turbulence 
that underlies the deviation angle statistics in the 
traditional interpretation of interstellar scintilla- 
tions. In contrast with the traditional picture, the 
cumulative effect of many independent scatterings 
by different thin shells is not Gaussian (which can 
be understood formally from the fact that the cen- 
tral limit theorem does not apply to slowly declin- 
ing distributions); it is instead a Levy flight with 
(3 ~ 2/3, as implied by the observations (see § 2). 

The foregoing derivation of the statistical prop- 
erties of scattering by thin shells is based ex- 
clusively on geometric considerations. If, how- 
ever, the shells are turbulent, then the effects of 
turbulence need to be considered in conjunction 
with those of geometry. Turbulence in the ionized 
boundary of a molecular cloud can be a natural 
consequence of the turbulent structure of the cloud 
itself or of the surrounding ionized medium. It is, 
however, difficult to construct a detailed model of 
how the turbulence is excited in the cloud or pene- 
trates into it since the theory of compressible tur- 
bulence is still in its infancy. However, the spec- 
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trum of a homogeneous, compressible turbulence 
has been investigated in some detail and found to 
be close to the Kolmogorov one (Boldyrev 2002; 
Boldyrev et al. 2002a, b). This should suffice for 
our qualitative estimate. 

We denote the outer scale of such a preexisting 
turbulence by i? ou t ■ This value should not be con- 
fused with the outer scale of turbulence excited 
inside expanding HII regions, which we analyzed 
in § 3. We assume that the density fluctuations on 
this scale are 5n out ~ n ~ 10 2 cm~ 3 (e.g., Larson 
1981). To simplify the analysis, we neglect the in- 
terplay of turbulent fluctuations with the bound- 
aries of the ionized shells. We consider the con- 
tribution of N slabs of homogeneous turbulence, 
each of width AR, that lie along the line of sight. 
The integrated density fluctuations are given by an 
expression analogous to equation (4), AjV(y) ~ 
no(y/i?out) (1+Q)/2 i?out 1/2 (iVAi?) 1 /2. When the 
interstellar ionizing radiation is not very strong 
(J ~ Jo), the boundaries are not fully ionized 
and ion-neutral collisions are important. One can 
assume that turbulent fluctuations are smoothed 
out below the collisional cutoff scale A m f p ~ 1.5 x 
10 15 /no ~ 10 13 cm, which is assumed to be domi- 
nated by collisions with helium atoms (e.g., Lith- 
wick & Goldreich 2001). This estimate indicates 
that we can safely assume that the turbulence is 
predominantly in the inertial range (with a w 2/3) 
on the scales of interest to us. 

Observations have revealed that the scaling of 
pulse broadening with dispersion measure changes 
from t cx DM for distant pulsars tor a DM 2 for 
nearby ones. The transition between the two scal- 
ing behaviors occurs at DM ~ 30 pc cm~ 3 (see, 
e.g., Fig. 4 in Bhat et al. 2004) and appears as 
an "elbow" in the t vs. DM diagram. This dis- 
persion measure corresponds to a pulsar distance 
d ~ 1 kpc, which, by the arguments given in § 3, 
implies a density fluctuation scale y ~ 10 10 cm. 
We propose that this change in the scaling be- 
havior is a reflection of the fact that, for nearby 
pulsars, the scattering angle is small and the scat- 
tering is dominated by turbulence within the ion- 
ized cloud boundaries rather than by the shape of 
these boundaries. Based on the arguments pre- 
sented in § 2, the correlation scale y of transmit- 
ted waves due to turbulent fluctuations satisfies 
the condition A<f>(y) = \r AAf(y) ~ 2tt. Setting 
A = 30cm, n ~ 10 2 cm~ 3 , and y ~ 10 10 cm, we 



infer a value of i? ut ~ 100 pc for the outer scale of 
the turbulence. It is gratifying to note that this es- 
timate agrees with the value commonly attributed 
(on both theoretical and observational grounds) to 
the general interstellar medium (ISM) turbulence 
(e.g., Armstrong et al. 1995). 

In the case of a strong turbulence inside the 
boundaries (J » J ) it is conceivable that the 
properties of the turbulence itself may play a key 
role in establishing the non-Gaussian nature of the 
scattering for distant pulsars (the first scenario 
outlined in § 1). In § 5 we consider the possi- 
bility that molecular clouds that are situated in 
the vicinity of a bright star or near the Galactic 
center (so that J ^> Jo) may naturally develop tur- 
bulence that is strong and non-Gaussian on small 
scales. We caution, however, that the properties 
of the envisioned turbulence are not yet well es- 
tablished. Furthermore, to fully analyze this case, 
one needs to specify how the turbulence is gen- 
erated and how it interacts with the boundaries 
of the ionized layer. The answer to both of these 
questions may depend on the particular mecha- 
nism of cloud formation, which may not be uni- 
versal and which, in any case, is beyond the scope 
of the present discussion. Therefore, in contrast to 
the robust and general conclusions that we reached 
in this section on the basis of purely geometric ar- 
guments, the results we obtain in § 5 are at this 
point still tentative. 

5. The Onset of Turbulence in Strongly 
Ionized Cloud Boundaries 

We consider only the extreme case, when a 
bright star is turned on very close to (or, alterna- 
tively, inside) a molecular cloud. In this case the 
surrounding gas gets ionized up to the Stromgren 
radius, R s , which can be evaluated from the ion- 
ization balance condition S = |7ri? 3 n|/32 (where, 
again, «2 is the density of ionized gas, fii is the re- 
combination coefficient, and spherical symmetry is 
assumed for simplicity). During the initial phase 
of fast ionization the gas density does not change 
(n 2 ~ no) and one obtains R s ~ 2 x 10 20 n 2 ^ 3 cm 
for S = 10 49 s _1 , which for no ~ 10 2 cm' 3 yields 
R s <~ 10 19 cm. The duration of the fast ioniza- 
tion phase is ~ 10 3 — 10 4 yr from the turn-on 
time of the ionization source (e.g., Spitzer 1978; 
Dyson & Williams 1997). The pressure imbal- 
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Fig. 3. — Sketch of the structure of an unstable 
ionization front. The width of the turbulent re- 
gion is of the order of the recombination length in 
the ionized gas, Z roc . The velocity of the uncom- 
pressed neutral gas in the reference frame of the 
shock is denoted by Vo, and that of the ionized 
gas behind the shock by V 2 . For a D-critical front 
both velocities are of the order of the sound speed 
in the ionized gas, C 2 . 

ance between the gas inside the Stromgren sphere 
(which is rapidly heated to T ~ 8000 K) and 
the surrounding cold molecular gas (T ~ 30 K) 
then leads to an expansion of the ionized region, 
which proceeds at approximately the sonic speed 
C 2 ~ 10 6 cm/s. This "intermediate" (or sonic ex- 
pansion) phase lasts for ~ 10 6 — 10 7 yr, until pres- 
sure balance with the ambient gas is restored at 
the final phase of the H II region's evolution. 

As the lifetime of the ionizing star is typically 
~ 10 6 yr, only the intermediate phase is relevant 
in practice. This phase nominally starts when the 
ionizing flux that reaches the boundary of the HII 
region drops below the critical value 2noC 2 . Af- 
ter this time a hydrodynamic shock forms ahead 
of the ionization front (IF). This shock compresses 
the neutral gas to a density n\ such that a weak D- 
type IF is set up (e.g., Spitzer 1978). In this config- 
uration (shown schematically in Fig. 3) the speed 
of the IF with respect to the compressed gas is 
smaller than the D-critical speed Vd = C\ /(2C2), 
where C\ <~ 3 x 10 4 cm/s is the sound speed in the 
neutral gas. This is equivalent to the condition 
J/ni < Vd, which implies that the neutral gas is 
compressed to m > JC2/C1, which is > 10 5 cm~ 3 



for J > Jo. 

In the limit when the relative speed of the IF 
and the compressed gas equals Vd , the ionized gas 
flows out of the front at the sound speed C 2 ■ Since 
this is also roughly the speed at which the front 
propagates into the (stationary) ambient gas, one 
can argue that the front should always be close 
to the D-critical state (Kahn 1954). Qualitatively, 
when the ionizing flux reaching the front drops 
below its D-critical value, the density of the ion- 
ized gas behind the moving front drops too. Since 
the photon absorption is proportional to n|, this 
would lead to an increase in the flux reaching the 
IF, which would turn it back into a D-critical front. 
Modeling and numerical simulations confirm that 
the speed of the ionized gas leaving the IF is in- 
deed near sonic, or possibly smaller by a factor of 
a few (Dyson & Williams 1997; Williams 2002). 

Weak or critical D-type IFs are subject to an 
instability with a linear growth rate ~ kC\ (for 
perturbations with wavenumber k), so the short- 
est modes initially have the fastest growth (Kahn 
1958; Vandervoort 1962; Axford 1964; Williams 
2002). The range of scales of unstable pertur- 
bations in real fronts is limited at the lower end 
by the front's width ijf ~ l/(aoni), which is 
< 10 12 cm in the case we consider, and at the 
upper end by the recombination length I 

rcc ~ 

C2 / \P2n2) , which for our fiducial parameters is 
~ 5x 10 16 cm. Numerical simulations by Williams 
(2002) have confirmed that this instability indeed 
develops and evolves into the nonlinear regime. 
It was found that an unstable IF becomes rough, 
with the perturbation wavelength and the pertur- 
bation amplitude saturating below the recombina- 
tion length. Wc conjecture (following the original 
speculation by Axford 1964) that the IF instability 
would generate strong small-scale turbulence. One 
can, in fact, argue quite generally that, in anal- 
ogy with a flow through a grid, a large-Reynolds- 
number outflow from a rough surface should be 
turbulent. The numerical simulations of Williams 
(2002) have indeed revealed the existence of large 
velocity shears on small scales, but the resolution 
of these experiments was not high enough to study 
the details of the indicated fine-scale structures. 
We note, however, that a compressible flow con- 
verging on a scale I will in general create density 
fluctuations of the order An 2 /n 2 <~ (AV r /C 2 ) 2 /2 
on the same scale. For the near-sonic conditions 
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characterizing the downstream IF flow, the den- 
sity fluctuations could thus be close to (within a 
factor of a few of) n-2- A schematic picture of the 
envisioned turbulence is shown in Figure 3. 

We propose that the turbulence likely excited 
behind IFs that propagate into strongly ionized 
clouds could generate the small-scale density per- 
turbations inferred in the Levy model of interstel- 
lar scintillations. A strong ionizing flux gives rise 
to a D-type front that produces a significant com- 
pression of the upstream neutral gas, and this in 
turn leads to a sufficiently large value of the ra- 
tio l T ec/ki oc 7ii/n2 for the turbulence to develop. 
(As we estimated in § 4, the ratio of the thick- 
ness of the ionized skin of a molecular cloud to 
the mean free path of an ionizing photon is only 
< 10 in the case of the typical interstellar flux 
Jo.) The estimated outer scale of this turbulence, 
^out ~ 10 12 — 10 16 cm, is much smaller than the 
<~ 10 pc scale that characterizes the turbulent H II 
regions invoked in the standard interpretation of 
interstellar scintillations (see § 3). However, if (as 
can be plausibly expected) the magnitude of the 
velocity fluctuations on the outer scale is AV < C2 
also in this case then the perturbations should 
be much stronger on the scales ~ 10 8 — 10 10 cm 
of interest. Furthermore, in the present inter- 
pretation the turnover time of the largest eddies 
{tout ~ 10 6 — 10 10 s) is typically shorter than 
the cooling time t coo \, and since the ratio of the 
turnover time to the cooling time decreases with 
the eddy scale, the difficulty associated with the 
expected damping of the compressible and entropy 
modes when these two times become comparable 
(see § 1) would be avoided. 

To illustrate the above points, we estimate the 
expected density fluctuation on the small scale 
y g ~ 10 8 cm. The density fluctuations at the 
outer scale are Ari2 < 10 2 cm~ 3 . Assuming a Kol- 
mogorov scaling (as proposed in the Goldreich & 
Sridhar 1997 turbulent cascade model), Sn(y)/n ~ 
(y/lout) 1 ^ 3 ■ Adopting an outer scale of <~ 10 15 cm 
(for which radiative damping should be unimpor- 
tant), we infer 5n(y g ) < 0.1cm -3 . For compari- 
son, the value obtained by adopting l out ~ 10 pc 
is a factor of ~ 30 smaller, and, even more signifi- 
cantly, in the latter case the fluctuation might be 
strongly damped by radiative cooling. 

In the proposed picture turbulence is continu- 
ously regenerated in the layer (of width ^ rcc ) be- 



hind the IF. The energy comes from the ionizing 
radiation and is deposited into fluctuating streams 
of hot gas that leave the rough surface of the front. 
Without such regeneration, turbulence would de- 
cay on the relatively short crossing time l ICC /C2- 
It is also worth noting that the outer scale of this 
turbulence in not the effective Stromgren radius 
(i.e., the thickness of the ionized cloud boundary), 
as is sometimes assumed (e.g., Yusef-Zadeh et al. 
1994) but rather the recombination length in the 
ionized gas, which is much shorter. 

To account for the inferred non-Gaussian na- 
ture of the density fluctuations we make a con- 
nection with the fact that in the Levy-flight in- 
terpretation the fluctuations responsible for the 
observed radio-pulse broadening span a narrow 
range of scales, 10 8 — 10 10 cm (see § 3), which is 
close to the estimated scale of the proton gyrora- 
dius pi in the ionized regions of molecular clouds 
(see § 1). Now, the turbulent MHD cascade dis- 
cussed above corresponds to shear Alfvcn waves 
only on scales 3> p%. On scales that approach the 
ion gyroradius the Alfven waves become kinetic 
and develop an electric field component that is 
parallel to the magnetic field. The electrons, as 
the most mobile particles, rapidly respond to this 
field by moving along the magnetic field lines in 
the direction of increasing electric potential. This 
provides a mechanism for coupling density fluc- 
tuations in the plasma to the cascade of shear 
Alfvcn waves. One can show that the coupling 
is strongest for transverse wavenumbers k± satis- 
fying k±pi <~ 1 and that the Fourier amplitude of 
the induced density fluctuations may be estimated 
from Snk/n ~ e^k/ksT <~ (k±pi)5Bk/B , where 
4> is the fluctuating electric potential, T is the tem- 
perature, and ^ is Boltzmann's constant (e.g., 
Camargo et al. 1996). 

The above coupling mechanism provides yet an- 
other means of accounting for the inferred pres- 
ence of strong density fluctuations on small scales. 
Even more significant, however, is the fact that, as 
the ion gyroradius scale is approached, the turbu- 
lence itself gets modified. This happens because in 
the presence of relatively strong electron density 
fluctuations and of the corresponding perpendicu- 
lar electric field (see the above estimate of 5rik/n), 
the ion polarization drift velocity becomes compa- 
rable to the velocity of fluid fluctuations, and so 
ions can drift compressively across the magnetic 
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field lines. This leads to the coupling of density 
fluctuations to the Alfvenic fluctuations; a similar 
effect is at work in the so-called electromagnetic 
drift wave (or drift Alfven) turbulence, which has 
been investigated in the context of plasma fusion 
devices (e.g., Hasegawa & Mima 1978; Hasegawa 
& Wakatani 1983; Terry & Horton 1982; Hazel- 
tine 1983; Camargo et al. 1996; Krommes 2002). 
Interestingly, numerical simulations of such turbu- 
lence have revealed a rather high level of density 
fluctuations (in equipartition with kinetic fluctu- 
ations) and statistical properties that are clearly 
non-Gaussian (e.g., Terry et al. 2001; Graddock et 
al. 1991). 1 

On the basis of the foregoing results one could 
plausibly attribute the generation of strong, non- 
Gaussian density fluctuations in turbulent, ion- 
ized molecular gas to the onset of compressive 
effects in a shear Alfven wave cascade that ap- 
proaches the proton gyroradius scale. (As we dis- 
cussed at the beginning of this section, IF insta- 
bility at the boundaries of molecular clouds pro- 
vides a likely mechanism for generating such tur- 
bulence that might by itself produce strong fluctu- 
ations on small scales.) More definitive statements 
must, however, await further studies of this turbu- 
lence in an explicitly astrophysical context. In the 
meantime one could perhaps look for a signature of 
the compressive effects in the density spectrum in- 
ferred from the pulsar scintillations data on scales 
<~ 10 8 — 10 10 cm. A likely signature would be a 
flattening of the spectrum on these scales: a flat- 
tening of this type was observed in the electron 
density spectrum of the solar wind (Coles & Har- 
mon 1989), where it was, in fact, attributed to 
kinetic Alfven wave effects (Hollweg 1999). 

6. Conclusion 

To conclude our discussion, we recapitulate the 
main ideas and results of this work. We have 
been motivated by a recent comparison of radio 
signals from distant pulsars with an analytical 
model (Boldyrcv & Gwinn 2003a), from which 
a non-Gaussian power-law distribution of wave- 
scattering angles was inferred (corresponding to 
Levy statistics). It was also deduced that the 
density fluctuations responsible for this scatter- 

We are grateful to Paul Terry for drawing our attention to 
these numerical results. 



ing should be rather strong on small scales (~ 
10 8 — 10 10 cm; see § 3). These conclusions dif- 
fer from the predictions of the classical theory of 
scintillations, which are based on Gaussian statis- 
tics. 

The inferred non-Gaussian behavior of the 
scattering process indicates that the scattering 
medium has a different structure than the one 
envisioned in the standard model, where the an- 
gular deviations are attributed to the cumulative 
effect of a large number of weak-scattering events 
along the line of sight. In contrast, we propose 
in this paper that the observed scintillations are 
dominated by rare scattering events in spatially 
intermittent density structures. We identify these 
structures with thin and curved layers of ionized 
gas, which could plausibly be the ionized bound- 
aries of molecular clouds. In this picture, a ray 
propagating from a distant pulsar can experience 
significant scattering only when its trajectory in- 
tersects such a boundary. The angular deviation 
induced during such an encounter can be rather 
large, so a single scattering event may dominate 
the total angular deflection. The trajectory of the 
ray therefore has the character of a Levy flight. 

For a simple calculation, we modeled the ion- 
ized cloud boundaries as thin spherical shells. We 
found that the scattering-angle probability density 
function predicted by this model has the asymp- 
totic form P(A8) cx (A6)~ 5 / 3 for large angular 
deviations (see § 4). Remarkably, this is precisely 
the form predicted by Boldyrev & Gwinn (2003a) 
on the basis of a comparison of the Levy model 
with observations. Our analysis therefore strongly 
suggests that the scattering structures in the in- 
terstellar medium are likely to possess a shell-like 
morphology. 

The spatially intermittent nature of the electron- 
density distribution is an important ingredient of 
our theory. So far, analytical and numerical inves- 
tigations of interstellar turbulence were mostly re- 
stricted to homogeneous settings with periodic (or 
otherwise simplified) boundary conditions. How- 
ever, the actual interstellar electron density distri- 
bution is not homogeneous, and geometric bound- 
aries of scattering regions can play a dominant 
role in the line-of-sight integrated density fluctu- 
ations. Such boundary effects have been observed 
in aero-optical experiments, in which an optical 
beam propagated through a confined turbulent 
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region (e.g., Dimotakis et al. 2001). 

The distant-pulsar relationship t oc DM 4 be- 
tween the pulse time delay and the pulsar disper- 
sion measure, for which the above model provides 
a natural explanation, is observed to change in 
nearby pulsars to r oc DM 2 , the expected scal- 
ing for Gaussian turbulence. We have interpreted 
this change as a consequence of the fact that, for 
nearby pulsars, the scattering is dominated by tur- 
bulence within the ionized cloud boundaries rather 
than by the shape of these boundaries. We have 
found that the measured value of the transition 
DM (~ 30 pc cm -3 ) is consistent with the inferred 
value (> 10 20 cm) of the outer scale of the general 
ISM turbulence. 

Although, as we just summarized, the non- 
Gaussian, Levy-type statistics of scattering angles 
may have a purely geometric origin, we also dis- 
cussed (see § 5) the situation in which the scat- 
tering regions are strongly turbulent, so that the 
contribution from turbulent density fluctuations 
dominates the scattering. This may happen when 
a molecular cloud boundary is ionized by radiation 
from nearby stars or by strong interstellar radia- 
tion in the vicinity of the Galactic center. In this 
case a near-critical D-type ionization front devel- 
ops and propagates into the cloud. 

Ionization fronts of this type are known to 
be linearly unstable. In its nonlinear phase, 
this instability could lead to a strong turbu- 
lence on scales below the recombination length 
^roc ~ 10 16 cm and above the ionization front 
thickness kf ~ 10 12 cm (where the numerical val- 
ues correspond to a cloud of density ~ 10 2 cm -3 ). 
The outer scale of this turbulence is ~ Z rec , so 
the density fluctuations could be strong enough 
on small scales to produce the inferred scattering. 
We contrast this situation with the standard the- 
ory, where the outer scale of turbulence is of the 
order of the Stromgren radius of an H II region 
(R s <~ 10 19 cm). In this case the density is pas- 
sively advected by the Alfvenic cascade from the 
large outer scale to the small dissipative one. This 
cascade has to pass through the radiative cooling 
scale (~ 10 14 cm), where the density fluctuations 
could be significantly attenuated, so by the time 
it reaches the much smaller scales implicated in 
the scattering process the turbulence might be too 
weak to account for the observed scintillations. In 
the case of a turbulent ionization front the outer 



scale of turbulence is smaller than the cooling scale 
and this problem is avoided. 

We also noted that the small scales implied 
by the non-Gaussian interpretation happen to be 
close to the ion gyroscale in an ionized molecu- 
lar cloud. The turbulent cascade would become 
compressible on these scales because of kinetic ef- 
fects, with the density fluctuations becoming cou- 
pled to the Alfvenic ones. This is expected to 
produce a flattening of the density spectrum in 
this range, providing an alternative means of ac- 
counting for the relative strength of the fluctua- 
tions on small scales. Even more tantalizingly for 
our contemplated application, there have been in- 
dications from numerical simulations for a clearly 
non-Gaussian behavior of the fluctuations in this 
case. However, as neither the nonlinear develop- 
ment of the ionization front instability nor the be- 
havior of kinetic Alfvcn turbulence in astrophys- 
ical contexts have yet been studied in detail, our 
discussion of the properties of turbulent scattering 
layers must be regarded as tentative at this stage. 
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